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Abstract 

For a dynamical system we will construct various invariant sets starting from its conserved 
quantities. We will give conditions under which certain solutions of a nonlinear system are also 
solutions for a simpler dynamical system, for example when they are solutions for a linear dynamical 
system. We will apply these results to the example of Toda lattice. 
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1 Introduction. 

A particle moving in the Newtonian gravitational field is known as the Kepler problem. The equations 
of motion are 



II* II 3 ' 

where x £ K. 3 \ {0} is the position vector. 

It is well known that the motions for the Kepler problem are planar. If we consider the plane OX1X2, 
the equations of motion becomes 

x\ = 2/1 

Xi = 2/2 
2/2 ^ 



These equations are Hamiltonian with the standard symplectic form on M 4 and the Hamiltonian 
function H = i (y? + 2/1) , l , ■ From Kepler's second law we have another conserved quantity given 

by A = xiy 2 - x 2 yi- 

For a > 0, consider the following conserved quantity K = H + -\ A. After a straightforward compu- 
tation we will obtain that V K — if and only if y\ = -\x 2 , and y 2 = — -ttXi and ||(xi,X2)|| = a 2 . Equiv- 
alently, the set {VA' = 0} is equal with the set {(xi,x 2 , 2/1,2/2) | (^1,^2) 1 (j/ij 2/2) = 0, and \\{xi,x 2 )\ \ = 
a 2 , and \\(yi,y 2 )\\ = ^ and sgn(yi) = sgn(x 2 )} which is invariant under the dynamics and is filled 
with solutions that represent uniform circular motions moving clockwise. These particular motions are 
also solutions for the linear Hamiltonian system 



Xi 


= -x 2 


x 2 


= Xi 


2/'i 


= -2/2 


2/2 


= 2/i i 



where the Hamiltonian function is A. 

Given the above analysis, we can raise at least two questions. Is it true that for a dynamical system 
that admits conserved quantities, the set of points where the gradients of these conserved quantities are 
zero, is an invariant set? When do two Hamiltonian systems have common solutions? In what follows, 
we will give an answer to these questions. 
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For the first question, the answer is positive and it is given in section two where we will also discuss 
various generalizations of this answer. In section three, we will present the conditions under which we 
can give an answer to the second question. In section four we will illustrate these results for the example 
of Toda lattice. Detailed computations arc presented in the Appendix. 

2 Invariant sets 

Let / : K n — > M. n be a C q , q > 1 function which generates the differential equation 

x = f(x). (2.1) 

We suppose that equation (|2.1[) admits a C q vectorial conserved quantity F : M™ — > K fe with k < n. 
We denote by Pi, P\. the components of F. For s £ {0, k} we introduce the sets: 

M F s) = {x £ M" I rankVF(x) = s} (2.2) 
where VF(r) is the Jacobian matrix 

VF(x) = . (2.3) 

For r <E {1, ...,(?} we introduce the sets: 

N[ r) = {x£M n \d a F i (x)=0, Vie{l,...,fc}, Vae{l,..,n} ! , VI < r} (2.4) 

where we note d a Fi — dx d Fi dx — if a = (a%, ai). 

Remark 2.1 We observe that 

M[ 0) =Nf 1) ={x£R n \-^(x) = 0, i£{l,...,k}, j£{l,..,n}}. (2.5) 

The set {MFA with s £ {0, k} is a partition of R". A critical point of F is a point in K" at which 
the rank of the matrix VF(ai) is less than the maximum rank. A critical value is the image under F of 
a critical point. The set of critical points of F is 

M c ir = ut 1 M ( ^ 3 . (2.6) 

Using Sard's Theorem (see [S]) we have that the set of critical values F(M F ) is of fc-dimensional measure 
zero providing that q > n — k + 1. 

Remark 2.2 We also have the obvious inclusions NK-s C NK_ C ... C N^y 

Theorem 2.3 The sets are invariant under the dynamics generated by the differential equation 



Proof. Because F is a conserved quantity, we have 

P(*t(x)) = F{x), 

where <I>t : R™ — > M. n is the flow generated by (|2.1I) . Differentiating, we have 

VF(*t(a))V$ t (!c)=VF(s). (2.7) 
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As V$t(a;) is an invertible matrix for any x G K™ which is not an equilibrium point for (|2.ip . (see 
1 ), we have that 

rank\7F($ t (x)) = rank\7F(x), 

which implies the stated result. ■ 

As a consequence we will obtain the following well known result which was applied for studying 
invariant sets of various mechanical systems, see for example [5J. 

Corollary 2.4 The set of critical points of F is an invariant set of the dynamics generated by the 
differential equation (|2.1[) . 



Theorem 2.5 The sets are invariant sets for the dynamics generated by the differential equation 

KB. 



Proof. Let i G {1, ...,k}, I G {l,...,q} and a = {a.\, ...,a{) G {l,...,n} . We will prove by mathe- 
matical induction that 

n 

d a Fi(x)= E {V®t) a Mx)-4V®t) ai p l (x)d F i ($ t (x)) + S i , a (t,x), (2.8) 

i8i,...,A=l 

where /3 = (0x, ...,/?;) and Si ta (t,x) is a sum with the property: "all the terms contain a factor of the 
form d~tFi(<f> t (x)) with \ j\ < \". 

We will prove this result by mathematical induction with respect to I. Componentwise the relation 
(|2.7[) implies our result for 1 = 1. 

Let a' — (ai, ai, G {1, ...,n} l+1 where a = (ai, aj) G {1, . Using the induction 

hypothesis we have: 

&*Fi(x) = 7 -^— ( J2 (^t) a Mx)...(^ t ) aiPl (x)d^F i (^)) + ^-(^(t,i)). 

By a straightforward computation we obtain 
d 



( (^t) a Mx)---(^t) ai ^(x)d Fi^t(x))) 



n r, 

= £ (V$ t )a 1 / Jl (a:)...(V$t)a l A(a:) s 0"F<(* t (aO)+ 

A,..,A=i Q,+1 

+ E « ((v$ ( ) aift (4..(v*«)« i ftW)3^ i ($ ( W) = 

n 

E (V$ t ) ctlft ( 2 ;)...(V$ t )a ! A(x)(V$ t ) Qi+1 A +1 (x)9' 3 > l (<i>t(x)) + 
/ 3 1 ,...,A,A +1 =l 

We will note that 

S i>a ,(i,z) = E 75 ((V-ft) ai p 1 (x)-(y^t) a ^(x))d p F l (^ t (x)) + - (S itCt (t,x)). 

All the terms of Si, a >(t, x) contain a factor of the form d 1 Fi(& t (x)) with | — ^ | < I + 1, which had to be 
proved. 
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Let /3 = (J3i, pi) £ {!,..., n} 1 and V$ 4 be the inverse matrix of V$t(x), where x is not an 
equilibrium point for (|2.ip . Consequently, we have 

n 

d0Fi($ t (x))= J2 (V^)^^)...^)^^)^^^)-^^*^)]. (2.9) 

ai,...,Q( = l 

Also, we will prove by mathematical induction that the sets are invariant under the dynamics 
generated by the differential equation (I2.1[) . For j = 1 we have iV^ = M£n, which is an invariant 
set (see Theorem 12.31) . We suppose that for all j < I the sets NE\ are invariant. Let x £ ^n+i) ^ e 
arbitrary chosen. Using (|2 .9[) for / + 1-order of derivation and the induction hypothesis, we deduce that 
$ t (x) £ Nn + i)7 Vt. Summing up, NF-^ are invariant sets for all j € {1, q}. ■ 



3 Finding solutions using simpler dynamics 

Let F, G : R™ -> R be C functions and the differential equations: 

x = f(x,VF(x)) (3.1) 

and 

x = f(x,VG(x)) (3.2) 

where / : R 2 ™ — >• R™ is a C" 3 vectorial function. 

For i £ 1", we denote with $f (x) and $f(x) the solutions of (|3.ip and (|3.2[) with the initial 
conditions 3>jf(a;) = a; and 3>q (a;) = a:. We will introduce the following set 

£i = {a: e R" | V-F(ar) = VG(x)}. 

Theorem 3.1 If F — G is a conserved quantity for i3.1\) and x € E±, then for all t we have 

<S>((x) = $?(x). 

Proof. For L = F — G we have the equality E\ = Mks . By Theorem (|2.3[) the set E\ is invariant 
under the dynamics of (|3.ip . For x £ E\ we have VF($f (a;)) = VG($f (a;)) for all i and consequently 

j t ^{x) = /(<Df (*), VF(*f (x)) = /($f (x), VG($f (a:)). 

The above equality shows that $f (a;) is also a solution for p.2[) . Given the uniqueness of the solutions, 
we obtain the desired equality. ■ 

Next we will discus a particular case of the result presented above. For this we will take f(x, y) — 
h(x) + g(y) and F = 0. Thus we have the two dynamics 

x = h{x) (3.3) 

and the perturbed dynamics 

x = h(x) +g(VG{x)). (3.4) 

We denote by $' t l the flux for (|3"3]l . 

Corollary 3.2 If G is a conserved quantity for (|3 . 3|) . then for all initial conditions x £ {x £ R™ | VG(x) = 
0} we have 
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Another particular case is when the function / in fl3.1[) and (|3.2[) verifies the equality 

< f(x,y),y >=0, 

where < •. • > is the Euclidian product on W 1 . In this case F is a conserved quantity for (|3.1[l and G is 
a conserved quantity for (|3.2j) . The above equality is verihed for the case when 



f(x,y) = U(x)y, 

where II(x) is an antisymmetric matrix. This is true for all almost Poisson manifolds (see [8j). In this 
situation the differential equations (|3.1I) and (13.21) become 

x = II(x)V-F(x) (3.5) 

and 

x = Il{x)VG(x). (3.6) 
We observe that F is a conserved quantity for (|3.6[) if and only if G is a conserved quantity for (|3.5j) . 



Corollary 3.3 If G is a conserved quantity for (|3.5| . then for the initial conditions in {x € K n | VF(i) 
VG(i)} we have $f (x) = $f (x). 



A particular case of Corollary 13.31 is when we have a symplectic manifold with G being a quadratic 
function. In this case, the solutions of the Hamiltonian vector field Xp starting in {x G K" | WF(x) = 
VG(i)} are also the solutions of the linear Hamiltonian system Xq. 

Analogous results are valid for the more general case of the vector valued conserved quantities and 
when the right side of equations (I3.1|) and (|3.2p depends on higher order derivatives. 

Firstly, we will introduce some notations. Let F, G : W 1 — > M. k be C q vectorial functions with q > 1 
and also k < n. If F\, F k are the components of F we denote with 

aIt-./ >, r°Fi , \ dFi dF 2 . . dF 2 . . dFk , . dFf. . k 
OTi aa; n ctei aa; n oxi aa; n 

and 

A r F(a;) = (...,a Q ^(a;),...) G R*"*, 

where r G {1, ...,«?}, i G {l,...,fc}, a G {l,...,n} r (|a| = r) and the components appear in the lexico- 
graphical order for (i,ce) in N r+1 . 

For a fix r G {1, q}, we will consider, as before, the two differential equations 

x = f r (x,A 1 F(x),...,A r F(x)) (3.7) 

and 

x = f r (x,A 1 G(x),...,A r G(x)). (3.8) 

Let x G R™, we denote with $f (x) and $p(x) the solutions of (|3. 7[) and (|3.8[) with initial conditions 
$q (x) = x and $g (x) = x. We will introduce the following set 

£ r = {x G K" | d a F(x) = d a G{x), V|a| < r}. 

Theorem 3.4 If F — G is a conserved quantity for (|3.7p and x G £v i/ien /or i we have 

<&f(x)=<I>p(x). 

Also the obvious extension of Corollary 13.21 takes place. 
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4 Invariant sets for Toda lattices 



The Toda lattice describes the one-dimensional motions of a chain of particles with nearest neighbor 
interactions. For a chain of particles with the equal masses m, Morikazu Toda came up with the choice 
of the interaction potential 

V(r) = e~ r + r- 1. 

The system of motion reads explicitly 

m ^ = e -(^-^-i)_ e -(^+i-^) j i€Z< (4.1) 

This second order differential system is equivalent with the first order differential system 

Xi = Ui , , 

where m is the velocity of the particle i. 

An equilibrium of the Toda lattice has the form 

Xi — + Xi, Ui = 0, A G M*, xo G R, i € Z. 

Let an equilibrium of the Toda lattice and = — x$ — Ai the displacement of the i particle from its 
equilibrium position. The system in the variables y, and Ui is 

Vi = u l x u on 

^ = £_( e -(j/i-»i-i) - e-Cw+i-fi)), i G z 1 ' J 



Let us define 



then the equations of motion become 



X, := e-GK+i-w), (4.4) 



Xj = Xi(ui -u i+ i) 

■di = Xi-x - Xi, ieZ. 



(4.5) 



The following particular cases are interesting: 

1. the case of a periodic lattice, X i+n — XiVi e Z, 

2. the case of a non-periodic lattice with the boundary conditions Xq — (correspond to formally 
setting yo = — 00 ) an d X„ = (correspond to formally setting y n +i = oo). 

In both cases we investigate the motions of the particles I to n (n 6 N* ) . 

4.1 The case of a periodic lattice 

In this case M. Henon proves in [5] that the following expressions are scalar conserved quantities 

Ira =J2 U ii" U id-Xh)-(-Xjl), (4-6) 

where m € {1, n} and the summation are extended to all terms which satisfy the following conditions: 

(i) the indices i%, ji + 1, + 1, which appear in the term (either explicitly, or implicitly 
through a factor Xj) are all different (modulo n); 

(ii) the number of these indices is m, i.e. k + 21 = m. Two terms differing only in the order of factors 
are not considered different, and therefore only one of them appears in the sum. 

In [4] , Flaschka has proved that the above functions are conserved quantities using a Lax formulation. 
This was generalized to arbitrary Lie algebras by Adler [2] and Kostant [7] . 

The first three scalar conserved quantities, depending on the variables (X%, ...,X n ,Ui, ...,u n ), are 

Ki<n 



G 



J 3 = u^u^u^ - ^2 UiXj, (Xq = X n ). (4.9) 

l<ii <i2<i3<ri j^i, j^i—l(modn) 

We will introduce the vectorial conserved quantities, 

Il2 = (Il,I 2 ), Il3 = (h,h), h 3 = (h,h), 1X23 = (hj2,h)- (4.10) 

I. The case rt odd 

For this case we obtain as invariant sets only subsets of the set of equilibrium points or the empty 
set. More precisely, Mfa = Alfa = Alfa = Mfa = Alfa = Alfa = Alfa 3 = Alfa 3 = 0. The following 
are subsets of the set of equilibrium points: 

Alfa ={(0,...,0,0,...,0)}, 

Mfa = {(x,...,x,o,...,o)\x €R}, 

Mfa = {(-^« 2 > -^u 2 ,u, «)|«e K}, 
Mg } 23 = {{X,...,X,u,...,u)\X,ueM}. 

II. The case n even 

We have, Mfa = Alfa = Alfa = Alfa = Alfa = Alfa 3 = Mfa 3 = Mfa 3 = 0. As nontrivial invariant 
sets we have the following: 

Mfa = {(Xx,X 2 , ...,X 1 ,X 2 ,ui,u 2 , ...,Ui,M 2 ) \Xi+ X 2 = uxu 2 , m + u 2 = 0}, 
Mfa = {(Xi,X 2 , Xi,X 2 , ui,u 2 , ...,m,u 2 ) I u\ + u 2 = 0}, 
Mfa = {{Xx,X 2 , ...,X 1 ,X 2 ,u 1 ,u 2 , ... , ui,u 2 ) \Xi+X 2 = -^(ui +u 2 ) 2 + uiu 2 }, 
Mfa 3 = {(X 1 ,X 2 ,...,X 1 ,X 2 ,u 1 ,u 2 ,...,u 1 ,u 2 )\X 1 ,X 2 ,u 1 ,u 2 e R}. 
We obtain Mfa 3 as the largest invariant set and the restricted dynamics is the dynamics of two particles 

r X 1 = X 1 {u 1 -u 2 ) 
I X 2 =X 2 (u 2 -u 1 ) 
I u 1 =X 2 -X 1 
{ ii 2 =X l -X 2 

On the invariant set Mfa we have the above dynamics subject to X\ + X 2 = — j(ui + u 2 ) 2 + uiu 2 . On 
the invariant set Mfa we have the above dynamics subject to u\ + u 2 = and on the invariant set Mfa 
we have the above dynamics subject to u\ + u 2 — and X\ + X 2 = u\u 2 . 

For the sets Mfa and Mfa the variables Xi have to take also negative values which from a mathe- 
matical point of view is correct and can be the solutions for the system ((4. 5|) . As the mechanical system 
is given by (|4.3|) and we do the change of variables (|4.4[) . the variables Xi have to be strictly positive 
in order to have a physical meaning. Consequently, from a mechanical point of view only the sets Mfa 

and Mfa 3 are meaningful. The computations can be found in the Appendix. 
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4.2 The case of non-periodic lattice 

It is known that if we have the matrices 



/ 




Xi 









( 


-Xi 









1 


u 2 X 2 












-x 2 ... 

















and B = 















1 


Mn-1 X n -i 












— X n - 


-1 


\ 





... 


1 u n 


) 




V o 








/ 



then the system (j4.5[) , in this case, has the Lax form 

t=[B,L], 



(4.11) 



where [B,L] = BL — LB. Using the Flaschka theorem (see [3]) we have the following scalar conserved 
quantities depending on the variables (X±, X ra _i, u\, ...,u n ) 



For k G {1, 2, 3} we have 



F k = -tr{L k ), fce{l,...,n}. 



Fi =J^Ui 

i=l 
n 2 

n— 1 ^ n 

f 3 = ^ Aj(uj + Ui+l) + o X/ 



(4.12) 



We will introduce the vectorial conserved quantities, 

F12 = (F 1 ,F 2 ), F 13 = (Fi.Fg), F 23 = (F 2 ,F 3 ), F 123 = {F 1 ,F 2 ,F 3 ). 
As before we will distinguish two cases. 



(4.13) 



I. The case n odd 

In this case we obtain as invariant sets only subsets of the set of equilibrium points or the empty set. 
More precisely, M% = M% = M?'f = M^f = M^ 3 = M^ 3 = M^ } 23 = M^f 3 =0. The following are 
subsets of the set of equilibrium points: 



M, 



F 3 _ 

(o) - 



M^ 3 ={(0,...,0,0,...,0)}, 



M: 



ffi = {(0, 0, ui, 0, ... Ul , 0, ui ) | ui 6 R} U {(0, 0, 0, u 2 , -, 0, u 2 , ) | u 2 € R}, 



II. The case n even 

We have, Af g = Afg - Af^ 2 
invariant sets we have the following: 



{(0, .■■,q, Mi,u 2 ^— I "1, "2 e 

n — 1 n 



A/; 



(i) 



(0) 



M: 



(0) 



Af; 



(0) 



M Fl23 
iw (i) 



As nontrivial 



Af ( ^ = {(A,0, A, 0, A, Mi, u 2 , ...,u 1 ,u 2 ) | mi + u 2 = 0, A = uiit 2 }, 
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Mf± 3 = {(X, 0, X, 0, X, ui, «a, ui.ua) I «i + u 2 = 0}, 

n— 1 n 

M&» = {(X, 0, X, 0, X, ui, u 2 , ui,u 2 ) | X = Ul u 2 }, 

W V v ' V * ' 

n— 1 n 

M'" s = {(X ) 0,...,X,0 ) X,u 1 ,u2,...,ui,u 2 )|X,ui,U2 € M}. 

1 ' v ' v ' 

n— 1 n 

We obtain M^ 23 as the largest invariant set and the restricted dynamics is given by 

X = Jf(u 1 -u 2 ) 
Ui = -X 
u 2 =X 

On the invariant set we have the above dynamics subject to X = uiu 2 . On the invariant set AffX 3 

we have the above dynamics subject to m + w 2 = and on the invariant set M^R we have the above 
dynamics subject to ui + u 2 = and X — uiu 2 . 

For the set Mj^? 3 the variables with i even are all equal with zero which from a mathematical 
point of view is correct. As before, the mechanical system is given by (|4.3p and we do the change of 
variables (|4.4[) . Consequently, the variables Xj have to be strictly positive in order to have a physical 
meaning. The computations can be found in the Appendix. 



5 Appendix 



The computations for the case of periodic lattice. 

We will make the notations 

l<i<n l<?i<?2<^ l^i^ n 

We observe that 

n 

2V = U 2 -J2^l (5-2) 

i=l 

With these notations we have the following, 

V/i = (0,...,0,l,...,l) (5.3) 

V/ 2 = (-1, -1, U_-ui, ...,U-u k , ...,U-u n ) (5.4) 
n+l n+k 2n 

V/ 3 = ~{U - u fe - u fe+1 ), V — u k (U — u k ) - (Y - X k _x - X k ), ...). (5.5) 

V v ' " v ' 

k n+k 

The study of M!* y 

The elements of M,^ are the solutions of the algebraic system 
U - Ufe - u k+ i = 0, 

V - u k (U ~ u k ) - (Y - X k _! - X k ) =0, Vfc g {1, n}. 



(5.6) 



Adding the first n equations, we obtain U — and u\ + u 2 = u 2 + U3 = ... = u n _i + u n = u n + u\. We 
deduce the following results: 
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The case n € 2N + 1. In this situation we have u% = u% = ... = u n = 0. Adding the last n equations 
of we obtain Y = and Xi + X 2 = X 2 + X 3 = ... = X„_i + X n = X n + X\. It implies that 

X\ = ... = X n = and consequently 

M£ } = {(0,...,0,0,...,0)}. 

The case n G 2N. For this situation Ui = (— l) l+1 u, u G K. In this case, using (|5.2j) . we have 
V = — § u 2 . The last n equations of (|5.6|) become 

-Y + AVi + X k = (~ - l) u 2 VA; G {1, n}. 

Making the addition we have F = — f w 2 and consequently X\ + X 2 = X 2 + A3 = ... = J„_i + A„ = 
A„ + Ai = — u 2 which implies Ai = A3 = ... = A„_i, and A2 = A4 = ... = X 2n . In this case we have 

Af(oj = {(Ai, A 2 , Ai, A 2 , iti,U2, ...,mi,M2) I Ai + X 2 = uiu 2 , wi + u 2 = 0}. 
The study ofMjj, M$ with G {(1, 2), (1, 3), (2, 3)}. 

A point (Xi, A„, Ui, u n ) G M?J? if and only if we have, for all fc, q G {1, n}, 



U -u k - u k+1 = 

F - u fe (C/ - u k ) - (Y - A fe _! — X k ) = V — u q (U - u q ) — (Y — X q _ x - X q ). 



(5.7) 



Adding the first n equations we obtain U — 0, u\ + u 2 = u 2 + 113 = ... = u n -i + u n = u n + u\ 
and consequently u% = {— l) l+1 u, u G M. The last n equations become Ai + X 2 = X 2 + A3 = ... = 
A„_i + X n = X n + X\ . 

The case n G 2N + 1. For this case note that u\ = U2 = ■■■ = u n = and Ai = ... = X n = A G K. In 
this case we have 

Mgj ={(A,...,A,0,...,0)|AgK}. 

The case n G 2N. It is easy to see that X\ = A3 = ... = A„_i, A 2 = A4 = ... = X n . In this case we 
have 

M {{) = {(Ai,A 2 ,...,Ai,A 2 ,ui,m 2 ,...,mi,m 2 )|ui +u 2 =0}. 
A point belongs to the set M?J? if and only if 

det(A feg ) = 0, det(B fcg )=0, det{C kq ) = 0, Vfc, g G {1, »}, (5.8) 



where 



B 



k(j 



A k = ( 1 1 

q \ -(U - U k - Ufc+l) —{U - U q - U q+ i) 

-1 U-U q 
-(U - life - u k+ i) V - u q (U - u q ) — (Y — A,_! - X q ) 



a 



kq 



U - Ufc U - u q 

V ~ u k (U - u k ) — (Y — Afc_x - A fc ) V - u q (U - u q ) — (Y — X q ^ - X q ) 



Using the expression of A kq we deduce that u\ + u 2 = u 2 + U3 = ... = + u n = u n + u%. 

The case n G 2N + 1. In this case we have u\ = u 2 = ... = u n = u, and U — nu, and V = ( ra ~ 1 )" M 2 
and 

-1 (n-l)u 



Bkq -\-{n-2)u ^-^-^ u 2 - {Y - X q - X q ^) 
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Using the condition that det(B kq ) = 0, wc obtain 

Y-X g - X q ^ = - -,r. v,/ m i » ! ■ 

Adding these relations we have 



(n-2)(n-l) , 



Y = 2~^ " 2 ' and X i = X 2 = = X n = X and X = -^-^-u 2 . 



With this relation we have 

n f (n — l)u (n — l)u 

kq = \ (n - 2)(n - l)w 2 (n - 2)(n - l)w 2 

We observe that the equality det(Ckq) = is verified. In conclusion we have 

^(i 2 ) = {(-^y^-w 2 , -^-u 2 ,u, u)|ue K}. 

The case n G 2N. In this case u\ — u 3 = ... = u n —x, u 2 = u 4 = .... = u n and we have U — + 1*2)- 
Using (|5.2|) we obtain V = "( n ~ 2 ) (u\ + u 2 ,) + ^u\u 2 . From the relation det(B kq ) — 0, we have 

n — 2 

y-.X,_i-.Y, = y-(l7-u,)(^— tf + u g ). (5.9) 

Consequently, Xi + X 2 = X 3 + X 4 = ... = X n _ 1 + X n and X 2 + X 3 = X 4 + X 5 = ... = X n + X x 
which further implies that X x = X 3 = ... = X n _ x , X 2 = X 4 = ... = X n and Y = + X 2 ). By 

substitution into (|5.9p we obtain X x + X 2 = — + u 2) 2 + uiu 2 . We have 

Afgj = {(X 1 ,X 2 , X 1 ,X 2 ,u 1 ,u 2 , ui,u 2 ) \Xi+X 2 = -f (mi + u 2 ) 2 + Ul u 2 }. 

The study ofMj™. 

We introduce, for k,q,r G {1, n}, the matrices 

1 

A kqr =| -1 -1 U — u r 

-(U - U k - ttfc+l) -([/- - Ug+l) U r 

1 1 

B kqr = I -1 U-U q U-U r 

-(U-Uk-Uk+l) V q V r 

/111 

Cfc ?r = U ~ U k U — u q U — u r 
V U fe U 9 U r 

where V, = V — u s (U — u s ) — (Y — X s _i — X,). We observe that ranfc(VIi23) = 2 <S> det(A kqr ) — 
det(Bkqr) = det(Ckqr) — Vfe, q,r £ {1, n}. The equations det{A kqr ) — Vfc, q, r G {1, n} give us 

iti + m 2 = u 2 + u 3 — ... = u„-i + it n = u n + u\. 

The case n G 2N + 1. We have u := u\ = ... = u n . With this notation we obtain U — nu, and 
V = ^^-u 2 and 14 = -^^-u 2 - (Y - X k ^ - X k ). It is easy to see that det(C kqr ) = 0. We have 
the equivalences 

1 

det(B kqr ) = det\ -1 (n ~ l)u ] = 

Ufc Vq-V r V r 
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V r — V q X r -\ + X r — X q -l + X q . 

Because n e 2N + 1, we obtain Xi — X 2 = ... = X„. In this case we have 

M^ 23 = {(X, X, u, m) \X,u£ R}. 
The case n e 2N. We have Ui = U3 = ... = 1*2 = M4 = ....= u n - By calculus we obtain 

/ 1 \ 

det(Ck qr ) = <=> <iet U — u k u k — u q u k — u r =0. 

V v k v q -v k v r -v k J 

If q = fc + 1 and r = fc + 2 we obtain Vk = Vk+2, which implies that X1+X2 = -X3+X4 = ... = X„_i+X„ 
and Jf 2 + X3 = X4 + X$ = ... = X n + X\. Consequently, we have X\ = X3 = ... = X n _\ and 
X 2 = X 4 = ... = X n . We observe that if s — t e 2Z, then u s — u t and V s = V tl which implies that 
det(Ck qr ) — Vfc, g,r£ {1, n}. For the matrices Bk qr we have the following properties: 

i) detBm = detBn2 = (by calculus). 

ii) detBkqr = —detB krq . 

iii) If qi — q 2 € 2N and ri — r 2 £ 2N then detB kiqiri = detB k2q2T2 . 

Using these properties we deduce that detB kqr = 0, Vfc, ij,re {1, n} and we obtain 

M^ 23 ={(X 1 ,X 2 ,...,X 1 ,X 2 ,u 1 ,u 2 ,...,u 1 ,u 2 )\X 1 ,X 2 ,u 1 ,u 2 e R}. 

The computations for the case of non-periodic lattice. 

If we consider the variables (X\, X n -\, u\, u n ) we obtain 

VFi = (0^,1^1) (5.10) 

n—1 n 

VF 2 = (1,...,1, «i,...,u„) (5.11) 

n-1 

VF 3 = {ui+u 2 ,...,u n - 1 +u n ,X 1 +ul,X 2 +X 1 +ul,...,X n - 2 +X n - 1 +ul_ 1 ,X n - 1 +ul). (5.12) 
The study of M^ 3 . 

F 3 
(0) 



The elements of are the solutions of the system 



f mi + u 2 = u 2 + U3 = ... = u n -i + u n = 

\ X 1 +u 2 1 =X 2 + X 1 + u\ = ... = X n _ 2 + X„_! + u\_ x = X„_! + u 2 n = 0. 
The case n e 2N + 1. Wc obtain 

M$ = {(0,...,0,0,...,0)}. 



(°) 

n—1 n 

The case n e 2N. In this situation ui = U3 = ... = it„-i = u, u 2 = 114 = ... = u n = —u and X k = —u 2 
if k e 2N + 1 and X k = if k e 2N. We have 

M (0) = {( ^, 0, 0,X , Ml,U2, . . .,Ml,U2 ) I "I + U 2 = 0, X = UlU 2 }. 

n — 1 n 

The study of M ( F «, M ( F ^ with e {(1, 2), (1, 3), (2, 3)} . 
The elements of are the solutions of the system 

( ui + u 2 = u 2 + u 3 = ... = u n -i + u n = 

\ Xi + u\ = X 2 + Xi + u\ = ... = X n _ 2 + X„_! + u\_ x = X„_i + u2. 
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The case n e 2N + 1. We obtain that 

M ( F - ={(0,...,0,0,...,0)}. 



n— 1 n 



The case n e 2N. In this situation 1*1=1*3 = ... = = u, 1*2 = U4 = ... = u n = — u and X^ = X 
if fc e 2N + 1 and X k = if k e 2N. We have 

Mjy = {( X,0,... , W,«b«2,---^b«2) I «1 + "2 = 0}. 
n—1 n 

For a point of the set M^™ we have 

cfei(A fc9 ) = 0, det(B fe<z ) = 0, det(C fcg ) = 0, Vfc, q (5.13) 

where 

* /V fe,«e{l,...,n-l} (5.14) 

u k + u fe+ i u q + u q+1 ) 



B * = {u k + u k+1 +X q+ u*)> - 1} and ? e {1, n} (5.15) 

X k - 1+ X k+ ul X q _ 1+ X q + ut)> *'«e{l, ...,»}. (5.16) 

Using the expression of we deduce that ui + it 2 = m 2 + m 3 = ... = u n -\ + u n and consequently 
u k = u q if k — q e 2Z. The matrices -B/c g have the form 

B kq =( ]_ „ "I , 2 ) if «62N+1 



and 



Bfe 9 = ( 1 y "I, , 2 ) if q G 2N. 

1 1*1 + M 2 Xg_l + X ? + M^ ' 



We have detB kq = if and only if Xi = Xi + X 2 = ... = X„_ 2 + -X1-1 = X„_i = U1U2. 

The case n e 2N + 1. It is easy to see that X\ = X2 = ... = X„_i = and M1M2 = 0. We observe that 

Mf™ = {(0, 0, mi, 0, u u 0, mi) I ui g K} U {(0, 0, 0, m 2 , 0, u 2 , 0) | m 2 e R}. 

\ > ✓ S ✓ ' ✓ N. ✓ 



n—1 n n—1 



The case n E 2N. We deduce that X l = X 3 = ... = X„_i = X, X 2 = X 4 = ... = X„_ 2 = and 
mim 2 = X. In this case we have 

M, F23 = {(X, 0, X, ...0, X, Ml, m 2 , Ml, m 2 ) I X = MlM 2 }. 

( ' S v ' V v ' 

n—1 n 

The study of M^ 23 . 

For a point of the set M^ 23 we have 

det{A kqr ) = 0, dct(B fe9r ) = 0, det(C kqr ) = 0, Vfc,<?,r (5.17) 

where 

1 

A kqr =\ 1 1 u r ] , k,q e {l,...,n- 1}, rg{l,..,n} (5.18) 

Mfe+M fe+ l Mq+Mq + l X r _l+X r +M,; 
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1 1 

1 u q u r 

Ufe + Ufe+i X q -i+ X q + u 2 q X r ^i+X r + ul 



, k e {l,...,n- 1}, q,r E {l,...,n} 

(5.19) 



1 1 1 

C kqr = ( u fe u 9 u r | , k,q,r e {1, ...,n} (5.20) 



Using the expression of A kqr we deduce that Ui + «2 = ^2 + «3 = ... = u n -% + u„ and consequently 
Uk = u q if fc — q £ 2Z. For fc G {1, n — 1}, we have 

/ 1 

detB k>k+ x >k = cfei 1 Ufe+i - u/c u/c 

\ Mfe + - Afc_i + W^+l - u\ X k -\ + X k + u\ 

and we deduce that X k +i = X k -i- 

The case n £ 2N + 1. It is easy to see that X\ — X 2 = ... = X n -% = 0, u\ =7/3 = ... = u n and 
U2 = U4 = ... = w„-i. All the conditions f|5 . 1 7(1 are verified and the set is 

M(2) 3 = {(p,...,0, Ml,M2 , -,Ml ) |«1,«2 e R}. 
n — 1 n 

The case n € 2N. It is easy to see that Xi = X3 = ... = X n _i = X, X2 = X4 = ... = X n -2 = 0, 
U\ — U3 = ... = and u 2 = = ... = it n . All the conditions (|5.17|) are verified and the set is 

M^f 3 = {( X, 0, ... , X, 0, X , ui,u 2 , u u u 2 ) I X, Mi, u 2 G K}. 

71 — 1 n 

Acknowledgments. Petre Birtea has been supported by CNCSIS UEFISCSU, project number 
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